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Abstract 

We present here a classical optics device based on an imaging architecture as analogy of a quan- 
tum system where the violation of the Bell inequality can be evidenced. In our case, the two qbits 
entangled state needed to obtain non classical correlations is encoded using an electromagnetic wave 
modulated in amplitude and phase. Computational states are represented in a way where each one 
of the two qbits is associated with two orthogonal directions in the input plane. In addition, unitary 
operations involved in the measurement of the observables are simulated with the use of a coher- 
ent optical processor The images obtained in the output of the process, contain all the information 
about the joint, marginal and conditional probabilities. By measuring the intensity distribution in the 
image plane we evaluate the mean values of the simulated observables. The obtained experimental 
results show, in an illustrative manner, how some correlations of Clauser-Horne-Shimony-Holt type 
exceed the upper bound imposed by the local realism hypothesis as a consequence of the joint effect 
of entanglement and two-particle interference. 

1 Introduction 

In this paper we present a classical analogy of some well known quantum experiments involving Bell 
inequalities violation. The analogy between quantum mechanics and classical optics has been recently 
explored [1-8]. The main idea is to exploit the wave nature of the electromagnetic field in order to rep- 
resent the quantum state of one or more particles. In this representation, the probability amplitude of 
ocurrence of each state of a basis is associated to the complex amplitude of the electromagnetic field and 
temporal evolutions are simulated by means of the propagation of the field through an optical system. 
Quantum phenomena can be understood as a consequence of the wave nature of the evolution of quantum 
states. In this sense, the wave character of the electromagnetic field allows us to simulate, in a pictorical 
way, the behaviour of the quantum world that usually contradicts common sense. Moreover, the classical 
electromagnetic field works as an ontologic representation of the wavefunction and it is a useful tool to 
visualize the structure of problems that ai^e usually complex and counterintuitive. 



Violation tests of Bell inequalities have become a fundamental tool for experimentally proving the pres- 
ence of entanglement correlations in quantum systems of general interest in ai^eas of quantum informa- 
tion, quantum computation and foundations of quantum mechanics. In a seminal paper [9] Einstein, 
Podolsky and Rosen (EPR) established their argument of the so-called local realism hypothesis. Accord- 
ing to it, if we accept that certain properties of a measured system are present prior to and independent 
of the observation, then quantum mechanics is not a complete theory of Nature. Almost twenty years 
later. Bell [10] showed that, for systems composed by two spin 1/2 particles, measurements of some cor- 
related quantities should yield different results in the quantum mechanical case to those expected if we 
accept the local-realism criterion of EPR. Many experiments confirmed the quantum predictions using 
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Bell-like systems as entangled photons in polarization degrees of freedom [11-18], entangled photons 
in transverse momentum degrees of freedom [19] and entangled atoms [20]. More recently, a novelty 
simulation of Bell inequalities violation using RMN techniques was reported [21] 

In this paper we will show that an analogy of Bell inequalities violation can be simulated using an optical 
architecture similar to those used in optical processing. In the Bell experiment, two subsystems shai^e an 
entangled state of two qbits. In our case the two qbits state is encoded by means of the optical modula- 
tion of an electromagnetic wave. The experimental set-up consists in an optical processor with a phase 
grating in its Fourier plane. The phase grating is represented in a spatial light modulator working in 
phase mostly mode. The output intensity distribution is registered by a CCD camera and is then analized 
with an imaging software. We show experimental results for a representation of a entangled state for 
which a Bell inequality of Clauser-Horne-Shimony-Holt (CHSH) type is violated [22]. We also discuss 
the illustrative interpretation of this kind of simulations. 

The paper is organized as follows: in Section 2 we give a brief review of the basic general concepts 
of Bell inequalities violation. In Section 3 we present some considerations about optical simulations of 
quantum information processing by means of optical architectures. We show how quantum states can be 
represented as images and universal quantum gates can be simulated as coherent optical processors. In 
section 4, we show how all these elements can be combined in order to obtain the optical setup represen- 
tation of two separated observers, each one performing measurements of local observables in the Hilbert 
space of their own qbit. In Section 5 we present experimental results in which a contradiction of the EPR 
local realism hypothesis can be evidenced. In Section 6 we provide some conclusions and discuss the 
implicancies of this simulation. 



2 Bell inequalities 

We will use in this paper the following notation: a state of the two-dimensional Hilbert space H2 (or 
the qbits space) is denoted as a complex linear- combination of the two states of the computational basis 
{|0) , 1 1)}. A separable state of the 2^ -dimensional space of a composed system of 2 qbits //f"^ is denoted 
as the product ® \'^{B))g, where = aj \0)j + Pj j = A, B; is the quantum state 

associated to the qbit A or B respectively. In what follows we briefly describe the CHSH approach 
to Bell type experiments [22]. Let us suppose that two spacelike separated observers (Alice and Bob) 
share an ensemble of entangled states 1^')^^^ = (|0)^ |0)^ + (X" /\/2 (subindices A and 
B denote the observers, or equivalently the qbits they will each measure). Let us consider two pairs of 
local physical observables: A and A' for Alice; and B and B' for Bob. As usually, we will define these 
observables as: 

A = a-aA, A' = a - aA, B = (3 ■ as, B' = P' ■ as (1) 

where 3, 3', /3, ai^e unit vectors and a a and gb are vectors whose components are the Pauli matrices 
operating on the local subspaces associated to Alice and Bob respectively. If Alice and Bob make a 
random choice of one observable of their pair and perform a simultaneous measurement, they have four 
possible non local combinations: 

A®B, A(g)B', A'(g)B, A'(S)B' (2) 

The measurements of the expected values of such nonlocal quantities can be carried out by performing 
local measurements. For instance, let us suppose that we want to know the expected value of A ® B. 
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After finishing a set of many experiments where Alice measures A or A' and Bob measures B or B', they 
perform the empirical sub-ensemble average of the product of their results {A(gi B) = 1/N J2f=i ^iBi 
where Alice has measured A and Bob B. The outcome of each local measurement is Ai,Bi = ±1 and N 
is the number of experiments where the Alice's choice was A and Bob's was B. It is of most importance 
to notice that Alice and Bob measure in coincidence so that the measurement which Alice performs does 
not disturbe the result of Bob's (or vice versa). They can do the same with the three remaining nonlocal 
observables A B', A' ® B and A' B'. In the end, Alice and Bob will have four independent 
measurement statistics of four randomly chosen sub-ensembles of the ensemble of singlet states. Alice 
and Bob can then calculate the expected value of the following nonlocal quantity: 

= A(g)B + A(S>B' + A'(g)B'-A'(g)B (3) 

Since PauU matrices do not commute, quantum mechanics asserts that the four nonlocal observables 
of eq.lEl are not compatible and therefore cannot be simultaneously measured. Moreover, the Bohr 
principle of complementarity [23] claims that we are forbidden to consider simultaneously the possible 
outcomes of mutually exclusive experiments. However, according to the local realism hypothesis, there 
exist local hidden parameters which completely determine the outcomes of the chosen measurement, 
having predictions quite different to those predicted by quantum mechanics. Moreover, these parameters 
also determine the outcomes that we would have if we have measured an observable which is incom- 
patible with that actually measured. We have no control on hidden parameters and so there are some 
degrees of freedom that are not precisely known. In the case of the Bell experiment, each hidden param- 
eter assigns well defined outcomes ±1 to each local measurement. It has been demonstrated that under 
such condition the expected value of the nonlocal quantity of Eq.[3] satisfies the CHSH form of the Bell 
inequality [22]: 

1(0)1 = \{A(S)B) + {A(g)B') + {A'(S)B')-{A'(S)B)\ < 2 (4) 

As mentioned in the introduction, quantum experiments measuring different observables have been per- 
formed [11-18]. We will inspect the particular case a' = P' = z so that A' = B' = C. With this 
assumption, {A' B') = (C ® C) = z • z = +1 for the state we have chosen, and the left Bell inequal- 
ity [4] becomes: 

(O) = {A® B) + {A®C) - {C ® B) <l (5) 

The equation above must be satisfied for arbitrarly chosen observables A and B with C = az- However, 
without loose of generality and for practical reasons, we will set observables whose optical implementa- 
tion is relatively simple as it will be clear in section 3. These observables aie: 

A = ax, B = smOcTx + cosOaz, C = az (6) 

Optical simulations of measurements of ax and az by means of imaging proccesing architectures have 
been implemented in previous works [6-8]. We will show in the next section that measurement of 
sin Oax + cosOaz can be performed in a similar way by controlling additional experimental parameters. 

Additionally, we will emphasize the rol of two-particle interference effects in the mechanism of Bell 
inequalities violation in this classical analogy. In order to make the rol of such interference effects to 
become apparent, we will analyze both the case of maximally entangled states and the case of mixed 
states. Let us consider the expectation value of the quantity O with respect to the mixed state whose 
density matrix is the convex sum of the pure entangled state Ppure = \^) ab {^\ab ^"^^ ^^'^ maximally 
mixed state: p^ixed = \ [(|0)^ |0)^) ((0|^ ® (Ol^) + (|1)^ ll)^) ((1|^ ® Therefore the 

state can be written as: 
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where the matrix representation is given in the computational basis {|0)^ , |1)^} {|0)^ , |1)^} and 
q G [0, 1]. The state of eq. [7] is a pure state for q = 1 and becomes mixed if < < 1 owing to the lost 
of the off-diagonal coherence components of the density matrix. The maximally mixed state corresponds 
to g = 0. Evaluating the expected value of the quantity O with respect to the state [7] and setting the 
observables defined in [6] we have: 



This last result contradicts the Bell inequality [8] if 6c{q) < < -k where 6c (q) is the solution of the 
equation q = {1 + cos 9) / sin 9 that always exists in the interval < < vr for all < g < 1. Maximal 
violation occurs for g = 1 and con^esponds to the Bell state which is pure and maximally entangled. 
On the contrary, there is no quantum correlations in the case of the maximally mixed state q = and 
therefore the Bell inequality is not violated. In the next section we will describe the problem underlined 
above with elements of classical optics. 



3 Optical simulation 

In order to perform the optical simulation of the Bell experiment, we will present first the way of encod- 
ing the complex amplitudes of the qbits as spatial distributions of light and we will emulate the unitary 
evolutions by means of an optical system composed of common optical devices such as lenses, phase 
gratings and phase shifters. These devices modify the complex amplitude of the electromagnetic field 
arising from the input image, where the state is encoded, reproducing the temporal evolution of the quan- 
tum state. The final state is obtained as an image that can be registered by a CCD. The spatial distribution 
of the intensity generated by the final complex amplitude can be interpreted as the probability distribu- 
tion of finding a given state. 



3.1 Encoding qbits in optical scenes 

We use a method which is extensively discussed in the Uterature [2, 3, 5] based on the representation 
of qbits as position C-bits and constitutes one possible scheme of the so called "unary representation of 
quantum systems". According to it, in the input scene we encode the logical values |0) and 1 1) of a single 
qbit in two slices located in the left and right halves of the full plane respectively. Let {xo, Vo) be the 
coordinates of the input plane. As the representation is one-dimensional, we have traslational symmetry 
with respect to the Uo axis, and we can consider only the Xq coordinate. According to wavefunction 
formalism, the optical analogy suggests the following notation: 



{0{9)) = tr [{A 5) p] + tr [{A C) p] - tr [{C 0B)p] = qsm9- cos 9 



(8) 




(9) 
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where Rect{x) is a unit rectangle function that takes the value 1 if |x| < 1/2 and in other case. Eq. [9] 
describes unit amplitude transmittance in a rectangle of width b centered in x = ±a (b <^ a) where the 
computational state |0) is associated to the left rectangle and the computational state |1) is associated to 
the right rectangle as it is shown in Fig. 1(a). To encode the more general state of two qbits, modulation 
of the complex amplitude of the field in both slices is needed as we can see in Fig. 1(b). The optical 
analogy of the quantum measurement process is very simple. In fact, if the state |(^) = a |0) + /3 |1) is 
represented, the relative intensities of the slices will be precisely |a|^ and |/?|^ and, after renormalization 
|a|2 + |/j|2 = 1^ they can be interpreted as the probabilities associated to the measurement outcomes 
+1 and —1 in a projective measurement in the computational basis. Making a step forward in this 
classical analogy, a quantum measurement is equivalent to keeping the left slice with probability \a\ or 
the right slice with probability Moreover, the expected value of az on the state represented, can be 
calculated by extracting information from the image represented in Fig. [l.b] as (cr^) = tr{az |<^) {(p\ ) = 
(|a|2-|/3p)/(|ap + |/?|2). 





Figure 1 : Optical representation of the single qbit state, (a) Optical representation of the states of the computa- 
tional basis, (b) The input scene associated to the optical single qbit. The state a |0) + /3 |1) is represented by the 
two single "left" or "right" slices where the constants a and (3 are the complex amplitudes of the electromagnetic 
field in each computational slice 

The scheme described above is easily generalizable to represent two or more qbits. In order to represent 
the computational states of a second qbit we will define another dicotomic regions within each one of the 
previously defined regions. As we can see in Fig.2, we will use the convention of splitting the plane in 
up and down regions to accomplish this. We denote A the qbit encoded in the up-down direction and B 
the qbit encoded in the left-right direction. In this case the representation of the two qbits basis will be a 
two dimensional extension of Eqs.[9]. 







|l>^®|0)s 





(a) 




Figure 2: Schematic picture of the representation of two qbits states by using optical scenes, (a) Spatial organiza- 
tion of the input plane in order to emulate two qbits states, (b) Optical representation of the |0)^ (g) \Q) g state, (c) 
Optical representation of the general pure two qbits state. Gray level scale correspond to different amplitudes and 
phase modulations of the classical wavefront 
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3.2 Pure and mixed states: the optical approach 

In the previous section we have discussed the optical representation and the classical interpretation of 
the measurement process in the case of pure quantum states. We will briefly introduce here how we can 
extend these concepts in the case of mixed states. A mixed state with density matrix p = i P« I V'j ) (v'i 1 5 
Pi > 0, J2iPi = 1' consists in a set of pure states H^Pi) ,i = 1,2, ....,n}, each appearing with its 
respective probability pi. The strategy we use to simulate and measure mixed states is as follows: first we 
represent an image in the input plane that simulate the member of the ensemble \ ipi), which is a pure state. 
Then we calculate the expected value of cr^, (cr^), as we have discussed in section 3.1, then we multiply 
the result by pi, and finally we sum over all the results. In this way, since tr{azp) = Yl,iPi'^'^{'^z Wi) 
the same values as having a statistical ensemble, will be obtained. It is worth mentioning that following 
this procedure another analogy can be achieved using a temporal succession of images with a duration 
Tj. If the total time of the experiment is T, then we select the times so that T-JT = pi. In this case, 
one has to integrate the output images in the time T of the experiment. This could be realized using an 
optical element capable of modulating the light field in the input plane in a dynamical way. For practical 
reasons, we used the first strategy throghout our work. 

3.3 The optical U(2) operator 

The key of our representation is the possibility of simulate the unitary change of basis that maps the 
complex amplitudes of the input state in the computational basis onto the complex amplitudes of the state 
in the basis of eigenstates of the observables that should be measured. This simulation is accomplished 
by using phase gratings as filters in an optical processor architecture. In what follows we denote the 
field amplitude as depending on a single relevant coordinate due to the unidimensional character of the 
optical simulation of local U (2) operators. The simulation of U(2) operators acting on single qbit states 
works as follows. The quantum state is encoded in an input scene located in the previous focal plane of a 
spherical lens of focal distance /. For collimated illumination, this lens allows to obtain on its back focal 
plane the Fourier transform of the input scene, which corresponds to a field distribution whose relevant 
coordinate i?, xp. The relationship between the spatial frequency variable fx and the position coordinate 
xp on the Fourier plane is fx = xp/^f, where A is the wavelength of the light field. In the Fourier plane, 
a spatial filter of complex transmittance H {fx) is placed. A second spherical lens of focal distance / 
is placed so that its previous focal plane lies in the Fourier plane. This lens allows to obtain the inverse 
Fourier transform of the product between the Fourier transform of the input transmitance and the function 
H {fx)- The field amplitude in the output plane is the convolution between the input complex amplitude 
and the so called impulse response of the system that is defined as the inverse Fourier transform of the 
function H {fx) [24]. In our case, spatial filtering in the Fourier plane is performed by an almenary phase 
grating which is a squai^e wave phase modulation of amplitude < (/> < 27r and spatial period 2p. We 
denote the width of each square pulse as p and the position of the center of the pulse in the frequency 
domain fc- The complex transmittance of the filter is: 



where the function above is defined in fx G [— p, p] and is extended by periodicity for all fx G 3^?- The 
complex transmitance defined in Eq. [10] could be expanded in terms of pure harmonic components so 
that the function H{fx) is written as H{fx) = Z^ne z CnCxp [i^fxj- Taking into account only the 
three central diffracted orders (coefficients Co and C±i), the inverse Fourier transform of the equation 
above is: 




(10) 
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where we have defined the phase constants 7^ = | it ^/c. Let us suppose now that in the input scene we 
represent a complex linear combination of the two computational states {xo\ ^'m) = « {xo\ 0)+/3 {xo\ 1) 
according to Eq.[9]. The output signal {xi\ '^out) will be equal to {h * ^m)(a; ) corresponding to the con- 
volution between the input signal and the impulse response of the system defined in Eq. [11] evaluated 
in the coordinate Xi of the image plane. The result gives six terms that are the three principal diffracted 
orders of the two slices that represent the computational states. Under certain conditions that can be 
experimentally controlled, the expression can be simplified. In fact, if we choose a grating whose spatial 
frequency satisfies the relationship 2p = l/2a with respect to the spatial separation 2a of the slices, 
the separation of the diffracted orders in the final plane will be equal to the distance of the two shces, 
allowing the interference between them, as it is suggested in Fig. [3]. All diffracted orders located out of 
the computational regions are not registered. 

The input-output relation of the process in matrix form, using the identifications of Eq.[9], is expressed 
in the representation of the computational basis as: 




cos t - sin t;e^'^ 
-sin^e*"^ cos ^ 




(12) 



where (j) and 7 are real- valued. The process is schematized in Fig. [3]. The similarity between the 
general expresion of U(2) operators and the 2-parameter family of linear operators of Eq. [12] becomes 
evident. For instance, if we want to measure the observable a^, we must apply the Hadamard operator 
H [25], and then perform the measurement in the computational basis. This can be done by setting 
tan((/)/2) = 2/7r and fc = —p/2. The Hadamai^d-like operator ^/2Haz obtained in this way, can be 
transformed in the proper Hadamard operator by placing a phase plate Uz in the front of the first lens 
since af = 1 Vi after renormalization. Measurement of hermitian observables defined in the Hilbert 
space of the two level quantum system can be simulated in the same way. In the Table [1] we show the 
unitary change of basis associated to each observable of Eqs.[6], and the corresponding values of the 
parameters (j) and fc- A phase shift az in front of the first lens must be eventually included. 

Hermitian Unitary (j) fc phase 
Ox H 2Arctan {2/tt) -p/2 ct^ 
sm.ea^ + cosOgz e'^^^y/'^Gz 2Arctan {{2/TT)Tan {9 /2)) -p/2 Uz 
1 -_ -_ - 

Table 1: Hermitian observables, unitary change of basis and parameters of the optical simulation 
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Figure 3: Schematic picture of the optically simulated U(2) operation. Complex amplitudes a and f3 are mapping 
onto a' and (3' by means of a 4f coherent optical processor with an almenary phase grating in the Fourier plane. 

4 Optical implementation of the analogy 

The experiment reported in this section, works as a classical optics analogy of the Bell experiment. We 
will test our setup in two cases which correspond to the two states q = and q = 1 described in Eq.[7]: 

• Case I: q = 1. Alice and Bob will share an entangled pair. For simulating this, we encode the 
maximally entangled state |^)^^ = (jO)^® |0)^ + |1)^<8> ll)^) /\/2by uniformly illumination 
of the top-left and the down-right quarters of the full input plane (Fig 4.a). 

• Case 2: q = 0. The input state will be the statistical mixture represented by the density matrix 

® |0)b) ((0^ «> (0|b) + (|l)^ ® \1)b) UMa ® (Mb)] whose optical analogy 
is the uniform illumination of the top-left or the down-right quarters of the full input plane, each 
with probability 1/2. This is an incoherent superposition of the computational states and it should 
not have any quantum correlation (Fig 4.b). 



(a) (b) 




p=l p=l/2 p=l/2 

Figure 4: Optical representation of the maximally entangled state (|0)^ ® |0)^ + ® (a) and of 

the maximally mixed state i [(|0)^ (g) \0) ^) ((0|^ ® (Oj^) + (g) 11)^) ((1|^ ® (Ij^)] (b) as optical scenes. 

The complete optical set-up is schematized in Fig. [5]. An Argon laser source (A = 477nm) is filtered and 
then collimated with lens Lq. The coUimated beam impinges onto the binary mask Pi which represents 
the two qbit state that Alice and Bob use during the experiment. The input scene is placed in the previous 
focal plane of the lens Li (focal length 26 cm) that allows to obtain the Fourier transform of the input 
in its back focal plane, or Fourier plane. In the Fourier plane, we place the spatial filter for simulating 
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two local unitary operators. According to previous discussions, this can be done with the composition 
of two orthogonal almenary phase gratings as showed in Fig. [6]. Horizontal phase grating with phase 
modulation parameter (pA produce diffracted orders in the "up-down" direction in where Alice's qbit is 
encoded. Vertical phase grating with parameter (pB produce diffracted orders in the "left-right" direction 
associated to Bob's subsystem. The two dimensional almenary phase grating, whose phase modulation 
goes from to (pA + 4>b {mod27r), was programmed in a spatial light modulator (SLM). This device 
consists in a Sony liquid crystal display TV (LCTV) that combined with two polarizers (Pi and P2) and 
two quarter wave plates (QWPi and QWP2), acts as a mostly phase modulator [26]. The LCTV (model 
LCX012BL) was extracted from a commercial video-projector and is a VGA resolution panel (640 x 480 
pixels) with squai^e pixels of 34 fim size separated by a distance of 41.3 /um. The process is completed 
with the lens L2 (focal length 26 cm) , that allows to obtain the inverse Fourier transform. The final 
image in the output plane Pq is captured by a videocamera (CCD). 




Figure 5: Experimental setup for simulating Bell experiment as an imaging system. 

According to Table [1], an adittional phase plate (PP) that introduces a tt phase shift in the left-bottom 
and in the right-top quarters of the input wavefront must be included. Therefore, this phase plate does 
not affect the illuminated zone of the input plane and we can ignore it. Moreover, in the final image, we 
must take into account the inversion of the coordinates system whose senses are indicated with arrows 
on the Pi and Pq input and output planes respectively in Fig. [5]. The protocol of the full experiment is 
depicted in Fig. [6] and can be described as follows: 

Each local operation is simulated by using one of the two orthogonal almenary phase gratings. Horizontal 
modulation from to (pA simulate unitary operations on the Alice "up-down" encoded qbit while vertical 
modulation from OtO(f)B works equally for the Bob "left-right" encoded qbit. Alice measurement is fixed 
to ax (or Uz) and therefore, eventually she only applies a Hadamaid operator before projecting her qbit on 
the computational basis. Bob measurement is varying as sin 6ax + cosOaz for < ^ < 27r and the phase 
modulation (I)b{G) = 2Arctan {{2/7r)Tan {6/2)) (see Table [1]) has to be variable. In both cases the 
amplitude of the phase modulation is controlled directly by the spatial light modulator. Measurement of 
az is performed by orthogonal projection on the computational basis and no phase modulation is needed 
(the white squares in Fig. [7]). The output scenes corresponding to the three pairs of local measurements 
(one fixed C and two varying A ^ B and C B) are registered by the CCD and recorded for its 
posterior analysis. The analysis method and the corresponding results will be shown in the next section. 
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Figure 6: Detail of tlie protocol of the full experiment. 



5 Bell inequalities: experimental results 

Once the output scenes corresponding to eacii measurement are obtained, the mean value of the observ- 
able (O) of Eq.[8] can be easyly evaluated. In fact, the output image contains all the information about 
the statistical properties of the measurement for a two qbits system. Unitary changes of basis reduce the 
problem of the measurement of an arbitrary observable into the problem of measuring tr^ (8> o"^. This is a 
projective measurement in the computational basis {|0)y^ , 1 1)^^} iSD {|0)^ , The results of a projec- 

tive measurement in computational basis can be easy interpreted in terms of the distributed intensities of 
the output field obtained in our experiment. The underlying process in the analysis of the output images 
is depicted in Fig. [7]. In what follows P{Ai, Bj) with i, j = ±1 is the joint probability of, after a pro- 
jective measurement in computational basis, A = i and B = j; i.e. the measurement outcomes obtained 
by Alice and Bob were i and j respectively. According to elementary probability theory, this quantity 
can be evaluated as P{Ai,Bj) = P{Bj/Ai)P{Ai); where P{Bj/Ai) is the conditional probability of 
B = j knowing that the result A = i was obtained, and P{Ai) is the probability that the result A = i 
was obtained, independently of Bob. 

Let us suposse that Alice performs the measurement of her qbit. As we have discussed above, it means 
that she has to keep the up or the down half of the full plane with different probabilities. The statistical 
properties of Alice's measurement are described by the marginal probabilities deduced from the joint 
distribution integrated over the degree of freedom associated to Bob. In the formalism of quantum me- 
chanics, marginal distribution of Alice is defined by the reduced density matrix trsipAB)- From the 
optical point of view, partial trace over Bob's subsystem means that the accesible information available 
for Alice is related to the field distribution in the "up-down" direction, independently to the field distri- 
bution in the "left-right" direction in where the information available for Bob is encoded. So, the two 
possible post-measurement states of the full system with their respective probabilities will be: 

|0)yi ® (aoo |0)b + "01 /\J\oim? + |aoi|^ (13) 



2 2 2 2 2 2 

with probability ^(^4+) = (I Qoo I + |aoi| )/(l«oo| + |aoi| + |«io| + |aii| ), that means that Alice 
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Bob post- measurement 
(conditional to the result of Alice) 

♦ 




P{A_,B_) = P{B_I A_)P{A_) 



Figure 7: Measuring a two qbits system in computational basis from the output distribution intensities, 
measurement outcome was +1 and her pseudorandom choice was the top half of the scene; or: 

(g) (aio |0)^ + an |1)^) /\J\aio\^ + |aii|^ (14) 

2 2 2 2 2 2 

with probability = (|aio| + |aii| )/(|aoo| + |aoi| + |aio| + |aii| ), that means that the 

outcome was — 1 and her choice was the bootom half of the scene. 

Now is the turn of Bob. The post-measurement state of the full system after Bob measurement will 
be conditioned for the result previously obtained by Alice. The four possible post-measurement states 
are naturally the four computational states. For instance, the post-measurement state |0)^ ® |0)^ can 
be obtained with probability = |aoo|^ /(|ooo|^ + |«oi|^) with the previous knowledge that 

the post-measurement state obtained by Ahce was the one described in Eq.[13]. At this point the non- 
local aspects of the joint measurement become evident. In this last case, the full process is a joint 
projective measurement of the output state in computational basis in where the result A = +1, 5 = +1 is 
obtained. The corresponding post-measurement state will be |0) ^ ^ |0) o with probability ^(^4-, -B+) = 

2 2 2 2 2 

P(S+/74+)i-'(A_|.) = |aoo| /(|«oo| + |aoi| + |«io| + |aii| ). The remaining three computational 
states, can be equally obtained as post-measurements states with probabilities depending on the intensity 
distribution \amn\ ; m, n = 0, 1. Therefore, the expected value of cx^ (go-^ in terms of the output intensity 
distribution is: 

{a,® a,) = P{A+,B+) + P{A^,B^)-P{A+,B^)-P{A.,B+) 
_ |«oo|^ + - |«oi|^ - |qio|^ 

I |2 I |2 I |2 I |2 v-'--^/ 

|aoo| +|aoi| + |aio| +|aii| 
It should be pointed out that Alice and Bob measurements involve only local operations. So the result of 
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the previous analisys does not depend on the chronological order of the measurements since local opera- 
tors commutes with each other. This is important since Alice and Bob measure simultaneously according 
to the EPR locality hypothesis. 



(b) 



(c) 



(a) 




Figure 8: First column: some typical output images. Second column: 3D plotting of intensity profiles, (a) Input 
state, (b) Output state after measurement of a z (Si ax- (c) Output state after measurement of ax ® (Jx- 

Calculations involved in Eq.[15] are performed directly from the relative output intensities associated 
to the computational states in the output image. Typically, images and 3D intensity profiles such as 
shown in Fig. [8] are obtained. The experiment consists in performing a sampling on the phase mod- 
ulation of the vertical grating with < (I)b < 27r, and for each value of evaluate the quantity 
(O) = {Ai^ B) + {Ai^ C) - {C B) in function of 6* = 2Arctan {{-K/2)tan (0b/2)) by inspecting 
the output images with the observables A, B and C defined as in Eq. [9]. The experimental results 
for (O) vs. 9 we. plotted toghether with the theoretical expected result gsin^ — cos 6. Theoretical 
and experimental curves are compared with 1 in order to explore possible violations of the inequality 
q sin 6 — cos 6 <1. Experimental results are summarized in Fig. [9]. 

In Fig. [9. a] (O) vs. 9 is plotted in the full range 9 G [0, 27r) for the maximally entangled state q = 1 and 
for the mixed state g = 0. As we can appreciate, although the experimental points differ slightly from 
the theoretical curves, a good cualitative agreement is obtained. It has to be mentioned that the SLM was 
used in maximal resolution of two pixels per spatial period of the grating. Since the properties of phase 
modulation of the SLM in high resolution are far from optimal, some significatively differences between 
experimental points and the theoretical curve appear, mainly in the range 9 £ [vr, 2tt) of Fig.[9.a]. Such 
differences ai^e more significant in the case q = 1 m where they are amplified by interference effects. 
Maximal violation of the Bell inequality occurs for q = \ and there is no violation for (7 = as expected. 
In Fig.[9.b] experimental and theoretical results corresponding to mixed states with density matrix de- 
fined in Eq.[7] are showed for g = 1, g = 2/3, g = 1/3 and g = 0. In this case the plot is shown not in 
the full range but in the zone of violation of the Bell inequality. The experimental results of the simula- 
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Figure 9: Experimental results, (a) Theoretical and experimental (O) vs. plotted in the full range 6 G [0, 2ti) 
for the maximally entangled state q = 1 and for the mixed state g = 0. (b) Theoretical and experimental results 
corresponding to mixed states with density matrix defined in Eq.[7] for q = 1, q = 2/?>, q = \/2> and g = in the 
Bell inequality violation domain. 



tion are also in good agreement with the theory, within the experimental error. The lost of coherence of 
the state, when the parameter goes from the maximally entangled to the maximally mixed state becomes 
evident in the transition from maximal violation to no violation. 



6 Conclusions 

We have implemented an optical setup to classically simulate Bell inequalities violation. We have shown 
how a conventional optical processing architecture can be used to optically simulate a Bell-type experi- 
ment scenario in where two space-like separated observers evaluate certain statistical correlations of non 
local observables by means of tensor products of local measurements. We have optically simulated all 
the elements corresponding to the real quantum process showing that all these elements have its classical 
optic counterpart. We have obtained experimental results by using the optical proposed set up. The sim- 
ulation begins with the optical representation of the quantum state of a four-dimensional Hilbert space as 
an image organized in four quarters according to the statement in Sec. 3.1. In this unary representation 
and according to Sec. 3.2, pure and mixed states can be emulated. The encoded information is orga- 
nized as a bipartite two level system. One part, conventionally called "Alice" has available dicotomic 
information related to the field distribution in the "up-down" direction of the scene. The second, "Bob" 
is associated to the ortogonal "left-right" direction. Quantum non locality arises from the assumption 
that the information available for one of the observers is unavailable for the other and vice versa. Then, 
we process the input image with an optical set-up composed by two sperical lenses and a bi-dimensional 
phase grating between them. The proposed architecture is in essential a 4f coherent optical processor 
with a phase grating in the Fourier plane. The optical processor is designed in order to ensure that the 
complex amplitude of the electromagnetic field is modified from the input to the output scene simulating 
an unitary evolution of the state. Simulated unitary evolution allows to measure tensor products of local 
observables. Mean values of such quantities are experimentally evaluated from the intensity distribution 
of the field in the output image. We show that, depending on the encoded input state, correlated quantities 
calculated from the expected values of the observables violate a Clauser-Horne-Shimony-Holt Bell-type 
innequality. 
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In order to emphazise the role of interference effects in the Bell inequalities violation, we test our set-up 
in two different cases: In the first, the maximally entangled state (|0)^ (8> |0)^ + \ ^)b) /\/2 en- 
coded in the input scene shows maximal violation meanwhile in the second, an incoherent superposition 
of |0)^ (g) |0)^ or (g) does not violate the Bell innequality mainly do to the absence of inter- 
ference effects. Parametrical lost of coherence from maximally entangled to maximally mixed states has 
been simulated by means of an optical representation of a convex mixing of both type of states. In all 
cases, the experimental results of the simulations are in good agreement with the theoretical predictions 
based on quantum mechanics. This simulation clearly illustrates how the mechanism of Bell inequalities 
violation needs both resources: quantum entanglement and two-particle interference. Even though the 
subject of this paper constitutes a classical simulation of the quantum reality that only has an illustrative 
interpretation, it nicely demonstrates that classical optics is a useful tool for a deeper understanding of 
some fundamental aspects of quantum mechanics sometimes considered to be hard and counterintuitive. 
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